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Abstract 

We generalise the definition of the characteristic of an integral triangle to integral sim- 
plices and prove that each simplex in an integral point set has the same characteristic. This 
theorem is used for an efficient construction algorithm for integral point sets. Using this 
algorithm we are able to provide new exact values for the minimum diameter of integral 
point sets. 



Key words: integral distances, minimum diameter 
2000 MSC: 52C10*, 11D99, 53C65 



1 Introduction 



Since the time of the Pythagoreans, mathematicians have considered geometrical 
objects with integral sides. Here we study sets of points in the Euclidean space 
E m where the pairwise distances are integers. Although there is a long history for 
integral point sets, very little is known about integral point sets for dimension m > 
3, see [3] for an overview. 

Due to Heron the area of a triangle with side lengths a, b, and c is given by 



(a + b + c) (a + b — c) (a — b + c)(—a + b + c) 

Aa 

4 

Thus we can write the area as = qy/k with a rational number q and a squarefree 
integer k. If ^ the integer k is unique and is called the characteristic or 
the index of the triangle. This invariant receives its relevance from the following 
theorem [4]. 
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Theorem 1 The triangles spanned by each three non collinear points in a plane 
integral point set have the same characteristic. 

This theorem can be utilised to develop an efficient algorithm for the generation of 
plane integral point sets [5,6]. Here we will generalise the definition of the charac- 
teristic of an integral triangle to integral simplices and prove an analogue to The- 
orem 1. Later on we will use this theorem to develop a generation algorithm for 
integral point sets in E m and present some new numerical data. 



2 Characteristic of integral simplices 

As the definition of the characteristic of an integral triangle depends on the area of 
a triangle we consider the volume of an m-dimensional simplex for point sets in 
E m . Therefore we need the Cayley-Menger matrix of a point set. 

Definition 2 IfV is a point set in E m with vertices vo,vi,..., i> n -i and C = (d% j) 
denotes the n x n matrix given by dfj = \\vi — Vj\\l the Cayley-Menger matrix C 
is obtained from C by bordering C with a top row (0, 1, 1, . . . , 1) and a left column 
(0,1,1,..., If. 

By CMD({v , v±, . . . , f„-i}) we denote the determinant of C({v , v±, . . . , v n -i}). 
If n — m + 1, the m-dimensional volume V m of V is given by 

K»W = ^de t (C). 

This allows us to define the characteristic of an m-dimensional integral simplex to 
be the squarefree integer k in V m (V) = qy/k whenever V m (V) ^ and q e Q. In 
order to prove the proposed theorem we consider a special coordinate representa- 
tion of integral simplices. 

Lemma 3 An integral m-dimensional simplex S = {v' , v[, . . . , v' m } with distance 
matrix D = (dij)o<ij< m and V m {S) ^ can be transformed via an isometry into 
the coordinates 



u = (0,0,...,0), 

V2 = (?2,l\/fcl, q2,2\fh, 0, . . . , 0) , 
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where ki is the squarefree part of v ^[^'^""'^ - )2 > %j £ Q, and qjj, kj ^ 

PROOF. We can obviously set v = (0,0,..., 0) and since e? ,i G N we can fur- 
thermore set vi = (do^y/ki, 0, 0, . . . , 0) where k\ = ^ferp = 1- Now we assume 
that we have already transformed v' ,v[, . . . , v\_ x into the stated coordinates. We 
set Vi = (xi, x 2 , ■ ■ ■ , x m ) with Xj E R. Since the points i> , ^i, . . . , i>j span an i- 
dimensional hyperplane of E m we can set — . . . — x m — 0. For j < % we 
have 

j , — i 

4i = H V i-^H2 = Y,^hh\J k h ~ X hf + X h- 
h=l h=j+l 

For < j < i we consider 



4i _ d h =J2 x l- (Qj,h\fk h - x h f 
h=i 

where we can set x h = q^h\fk~h for h < j by induction, yielding 

E 

Thus 

j'-i 



<t - 4* = _ ?J,A + 2( lj,3\fkh x 3 + 2 2 %hqj,hk h - ql h k h . 



Xj — 



Qjjh + E - 2qi,hqj,h k h) + d 2 0>i - dj. 

h=l 



Zqjjvhi 



and we can write = qijykj since 2qjj^/k h ^ due to induction. With this we 
have 

^o,j = E x h = x 1 + X 1i,h k h- 
h=l h=l 

Thus 



i-1 



h=l 



We also have q^^/kl ^ since i>q, . . . , v[ cannot lie in an i — 1 -dimensional 
hyperplane of E m due to V m (v' , v' 1: . . . , v' m ) ^ 0. □ 

The fcj are associated to the characteristic char(5) = k in the following way 
char(5) = k = squarefree part of JJ fc, . 

3=1 

Theorem 4 /n an m-dimensional integral point set V all simplices S — {v , v±, . . . , v m } 
with V m {S) 7^ have the same characteristic char(S) = k. 
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Proof. It suffices to prove that char(5i) = 01101(52) for two integral simplices 
Si = {vo,V!,....v m } and<S 2 = {%...,v m _i,i/ m } with V m (Si), V m (<S 2 ) 7^ 0. 
With the notations from Lemma 3 we have for the distance between v m and v' 



1=1 

m 

= 5Z(9m,jy ki — Q m i\J k{) + (c[m,m\l k m — (l m) m\j k' m ) 



i=l 
m—l 



i=l 



Thus yk m ,k' m has to be an integer. Because k m and k' m are squarefree integers 7^ 
we have k m = k' m and so char(5i) = char(5 2 )- □ 



3 Construction of integral point sets 



The key principle for a recursive construction of integral point set consisting of 
n points is the combination of two integral point sets V\ = {v , . . . , i>„- 2 } and 
V2 = {vq, • • • , t>n-3, v n -i} consisting of n — 1 points sharing n — 2 points, see 
Figure 1. Here we describe an integral point set by a symmetric matrix D = (d iy j) 




Figure 1. Combination of two integral point sets. 

representing the distances between the points. Because not all symmetric matrices 
are realizable as distance matrices in E m we need a generalisation of the triangle 
inequalities. 

Theorem 5 (Menger [9]) A set of vertices {vo, vi, . . . , i> n -i} with pairwise dis- 
tances d it j is realizable in the Euclidean space E m if and only if for all subsets 
{i , i r -i} C {0, 1, . . . , n — 1} of cardinality r < m + 1, 

(-l) r CMD({v io ,v il ,...,v ir _ 1 })>0, 

and for all subsets of cardinality m + 2 < r < n, 

(-iyCMD({v io ,v il ,...,v ir _ 1 }) = 0. 
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Fortunately we do not need to check all these equalities and inequalities. Because 
the point sets V\ and V 2 are realizable due to our construction strategy it suffices to 
check (-l)-CM%^,.,vi}) [5]. 

To solve the equivalence problem for integral point sets we use a variant of orderly 
generation [ 1 ,7,8, 1 1] . For the required ordering we consider the upper right triangle 
matrix of D leaving out the diagonal, 

^0,1 ^0,2 • • • do,n-l 
dl,2 ■ ■ ■ ^l,n-l 

1 

\ d ri _ 2 ,n-l J 

and read the entries column by column as a word 

w{D) = (do,i, d 0i 2, di,2, , do,n-i, ■ ■ ■ , d n _2, n -i) • 

With a lexicographical ordering on the words w(D) we define 

D 1 h D 2 wiDi) y w(D 2 ) 

for distance matrices Di, D 2 . We call a distance matrix D = (d i: j) <i : j <n canonical 
if 

D h (d T (i), T (j)) Vr e S n . 
By [D we denote the distance matrix consisting of the first n — 1 rows and columns 
of D. With this we call a distance matrix D semi-canonical if 

ID ^|(rf r (i),r(j)) VreS n . 

A canonical distance matrix is also semi-canonical. It is left to the reader to prove 
that each semi-canonical distance matrix D can be obtained by combining a canon- 
ical distance matrix D 1 and a semi-canonical distance matrix D 2 , see Figure 1. 
Only the distance rf„_ l n _ 2 is not determined by the distances of D x and D 2 . Here 
we consider two cases. If we combine two (m! — 1) -dimensional simplices to get 
an m'-dimensional simplex Theorem 5 yields a biquadratic inequality for G?„_ l n _ 2 - 
In the other case we can determine one or for n = m + 2 at most two different 
coordinate representations of the n points similar to the proof of Lemma 3, calcu- 
late (i n _i in _2, and check whether it is integral. We denote the sub routine doing this 
by combine(D 1 , D 2 ). At first we provide an algorithm to generate m-dimensional 
integral simplices. Therefore we assume that for a given diameter A, this is the 
largest distance, we have two lists C, c m , C s m of the canonical and the semi-canonical 
(m — 1) -dimensional integral simplices with diameter A which are ordered by -<, 
respectively. The following algorithm determines the lists C c m+1 and C s m+1 of the 
m-dimensional integral simplices with diameter A ordered by -<. 
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Algorithm 6 

Input: C c m , C s m 
Output: C c m+1 , C s m+1 
begin 

loop over x e £^ do 

loop over 9 y ^ ar with [x =[y do 
loop over z e combine(x, y) do 

if z is canonical then < — 2 end 
if 2; is semi-canonical then C s m+1 < — z end 
end 
end 
end 
end 



Because an m-dimensional simplex is an m-dimensional point set consisting of 
n = m + 1 points we can use Algorithm 6 to generate complete lists A4^ +1 , 
M-m+i of the canonical and semi-canonical m-dimensional integral point sets with 
diameter A consisting of m + 1 points, respectively. An m-dimensional point set is 
in semi-general position if no m + 1 points are situated on an (m — 1) -dimensional 
hyperplane. Using Theorem 4 we can give an algorithm to determine the lists M. c n 
and M. s n of the m-dimensional integral point sets in semi-general position consist- 
ing of n points with diameter A. 

Algorithm 7 

Input: M c n _ x , M s n _ x 
Output: M L n ,M s n 
begin 

M c n = ®, M s n = ® 
loop over x e M. c n _ x do 

loop over M s n _ x 3 y <x with [x =[y and char(x) = char(y) do 
loop over z e combine(x, y) do 

if z is canonical then M c n < — z end 
if z is semi-canonical then M s n < — z end 
end 
end 
end 
end 
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Table 1 



Number of calls of combine(x, y). 
4 Improvements 

To demonstrate the significance of Theorem 4 for an efficient enumeration algo- 
rithm for integral point sets we compare in Table 1 the number \l/(3, A) of calls 
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of combine(x, y) in Algorithm 7 for m = 3 and n = 5 to the number \I>(3, A) of 
calls of combine(x, y) without using Theorem 4. Additionally we give the number 
<5(3, A) of semi-canonical integral tetrahedrons with diameter A. 



5 Minimum diameters 



From the combinatorial point of view there is a natural interest in the minimum 
diameter d(m, n) of m-dimensional integral point sets consisting of n points. By 
d(m, n) we denote the minimum diameter of m-dimensional integral point sets in 
semi-general position. If additionally no m + 2 points lie on an m-dimensional 
sphere we denote the corresponding minimum diameter by d(m, n) and say the 
points are in general position. To check semi-general position we can use the Cayley- 
Menger matrix and test whether V m = or not. In the case of general position we 
have the following theorem. 



Theorem 8 Given m + 2 points in E m , with pairwise distances dij and no m + 1 
points in anm — 1-dimensional plane, lie on an m-dimensional sphere if and only 
if 





d 2 
u i,o 



d 2 

U 0,m+1 



d 2 

^m,m+ 1 



d 2 

"m+1,0 



d 



m+l,m 



See [2,10] for a proof. 



We have implemented Algorithm 6 and Algorithm 7 and received the following 
values for minimum diameters, see also [3,6,10]. The values not previously known 
in the literature are emphasised. 

d(3, n) 4 < n < 7 = d(3, n) 4 < n < 7 = 1, 3, 16, 44. 

rf(4,n) 5 <„< 8 = 1,4,11,14. 
d(4,n) 5 < n < 8 = 1,4,7,14. 
^(5,n) 6 <„< 9 = d(5,n) 6 < n < 9 = 1,4,5,8. 



To determine d(m, n) we have to modify Algorithm 7 because not every m + 1 
points of an m-dimensional pointset span an m-dimensional simplex. So we have 
to combine lower dimensional point sets with m-dimensional point sets. We leave 
the details to the reader and give only the results, 
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d(3, n) 4 <n< 23 = 1, 3, 4, 8, 13, 16, 17, 17, 17, 56, 65, 77, 
86, 99, 112, 133, 154, 195, 212, 228 
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